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ee ABSTRACT 


The stability of an infinite slab of an incompressible, inviscid fluid of low electrical con- 
ductivity carrying a uniform current produced by external sources is investigated in two 
eases: (1) when both the surfaces of the fluid slab are free and (2) when one surface is free, 
the other being rigid. In the latter case the influence of gravity acting perpendicular to the 
surface of the slab is included. The influence of surface tension of the fluid is, however, taken 
into account in both the cases. The effect of a homogeneous magnetic field produced by 
external sources and directed along the current in the fluid is also discussed. 

- It is shown that the fluid slab is always unstable when it carries an electric current which 
is closed outside the fluid symmetrically at infinity. In particular when the slab has both of 
its surfaces free the fundamental modes are either symmetric or anti-symmetric about the mid- 
plane in the absence of external magnetic field while they are neither symmetric nor anti- 
symmetric in its presence. 


1. Introduction 


Kruskal and Schwarzschild [1] studied among other problems the stability of 
a compressible fluid slab of infinite electrical conductivity with volume as well 
as surface current. The fluid is assumed to have the magnetic field directed 
everywhere perpendicular to the current and extending outside the fluid also. 
Assuming the strength of the field outside the fluid to be larger than the inside 
value, they find that the fluid slab is stable to sinuosoidal perturbations of the 
boundaries. A similar problem has also been studied by Loughhead [2], who 
showed that an,incompressible fluid slab can be unstable for some wavelengths 
if the magnetic field outside the fluid has a smaller value than inside. 

In the study of the stability of the fluid slab the above-mentioned authors 
assumed that the fluid had infinite electrical conductivity. This assumption is 
good in the case of a plasma and under the conditions one finds in astrophysi- 
cal problems. In these situations one rarely comes across a geometric configura- 
tion which can be approximated to a slab of finite thickness and extending to 


1 On deputation from the Tata Institute of Fundamental Research, Apollo Pier Road, 
Bombay 1, India. 
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ah i in the : nductivity bu 
d characteristic dimensions of the problem are of the order of 1 cm and there 
fore the assumption of infinite electrical conductivity is not justified. On th 
other hand, one can obtain the hydromagnetic stability criteria for fluids w 
electrical. conductivity is of the order that of mercury in the limit of zerc 
electrical conductivity if the currents are maintained by external surces [3, 4], 
It is the aim of this paper to study the stability of a fluid slab carrying uni 
form volume current in the presence of a homogeneous external magnetic field 
directed along the current under the assumption of low electrical conductivity 
of the fluid. ; e 

This problem will be studied in two parts. In the first part the stability of 
a fluid slab with one surface rigid and the other free, taking into account the 
effect of acceleration due to gravity acting perpendicular to the surface of the 
fluid, will be investigated. In the second part the stability of the fluid slab with 
both surfaces free, neglecting the effect of gravity, will be investigated. The 
influence of surface tension will be included in both the cases. 


ie 


PART I 


i i i ling  emet 


Instability of a fluid slab with one surface free and the other rigid 


2. Formulation of the problem 


Let us consider a slab of an incompressible, inviscid and non-magnetic fluid 
confined between two parallel planes separated by a distance 2d. Let us choose 
a right-handed Cartesian coordinate system with the z-axis perpendicular to the 
surface of the fluid with the origin situated midway between the planes. Then 
the two planes of the fluid slab are given by z= +d. The fluid is assumed to 
extend to infinity in x and y directions. Let there be an uniform current of 
constant volume density, j), flowing in the y direction inside the fluid slab and 
closed at infinity symmetrically on either side of the slab. The solution of the 
problem when the current is not symmetrically closed offers no difficulties, but 
we shall confine ourselves to the above assumption in the present paper. A 
uniform magnetic field, B,, produced by external sources is also present and is 
directed along the y-axis. We shall assume that the plane z= —d is rigid, while 
the plane z=d is free. The force of gravity is assumed to be acting in the 
negative z direction (see Fig. 1). We shall further assume that the regions z>d 
and z>—d are vacuums. 


The magnetic field distribution due to the current is given in rationalised 
M.K.S. units by 


B= (B,, 0, 0), (2.1 a) 
where B,=UMojoz if —d<z< +d, (2.1 b) 
By=pPojot, Wh 22d, (2.1 ¢) 
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so ait dies ie in Bcenbriins with eee iGiribution given inside 
gee 8 


p= koe *)— og (z—d), 


re Aa re of the fluid and g=acceleration due to gravity. 
order to study the stability of the fluid slab we assume that the free sur- 


ace ot the fluid is perturbed and its shape is then given by 
z=d+a cos (k,x+k,y), e 3) 


vhere ay The signs of k, and k, are for the present kept arbitrary. We 
ume further that a is periodic in time such that ; 


2 a=, cos wt. (2.4) 


The purpose of the theory is to obtain the dispersion relation between w and 
3 and k,, which enables us to find out the range of values of k, and k, for 
which @ has an imaginary value which corresponds to unstable modes of the 


perturbation. 


_ The basic equations of the problem are 

.* d 

y iV A a 

s : Cae wa Pe Best, (2.5) 

: VV" y=0, (2.6) 

% 

r Vx B= fyi, (2.7) 

= vy: B=0, (2.8) 

; I =constant. (2.9) 
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amperes. 


3. Solution for w _ a | 
(a) Current density distribution in the magnetostatic approximation 


As stated above, the current density is the gradient of a harmonic function. 
This fact together with eq. (2.9) and the boundary conditions that there should 
be no component of current density vector normal to the boundary given by 
eq. (2.3) and ‘also at the boundary z= —d enables one to determine’ uniquely 
‘the current density distribution inside the perturbed fluid slab. In the calcula- 
tion of any physical quantity we shall retain only terms proportional to a/d 
and neglect higher order terms. Following the method adopted in [8, 4] it can 
be shown that the perturbed current density is given by 


j= (Jes dus Je) (3.1) 

where jz =k, [Ay cosh kz + A, sinh kz] cos g, (3.2a 
Jy =o + ky [Ay cosh kz+ A, sinh kz] cos g, (3.2 b) 
and je=k[A, sinh kz+.A, cosh kz] sin g, (3.2.c) 


where A, and A, are constants given by 


jg ky 


os es Wiegand C= 

d Pig: jg ky 
eS Ay 2k cosh kd’ (4 
where k= + (k2 + 2), (3.5) 
and yp=(k,x+ kyy). (3.6) 


(b) Magnetic field distribution 


The distribution of the magnetic field inside the fluid slab is determined by 
eqs. (2.7) and (2.8). In the region outside the fluid, the ‘magnetic field vector 
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Biss : ow aes 4. Bil sane ¢ . j ra 
B 408 (4, sinh bet Ay cosh Ka)} 00s (3.7) 
By = s (C, sinh kz + C, cosh kz) 
Moke 4 «3 , 
as (A, sinh kz+ A, cosh koh Cos @, (3.7 c) 
pat : » 
and Bz=(C, cosh kz+ C, sinh kz) sin 9, (3.7 d) 


vhere C, and C, are constants given by 


1 : 
C,=C,= — 95 Moo ke exp (—kd). (3.8) 


e . (c) The pressure distribution 


Tt has been shown by Murty [3, 4] that the pressure distribution inside the 
‘uid can be obtained from the solution of a Poisson’s equation with the right- 
hand side given by the current density alone when the current density vector 
is a gradient of a scalar function and if there are no other forces acting on the 
fluid. However, the presence of a volume force due to gravity does not intro- 
duce any difficulty if the sources of the force are outside the fluid (i.e. if we 
‘neglect the self-gravity of the fluid). This can be seen as follows. If we take 
‘the divergence of both sides of eq. (2.5) and use eq. (2.6), we obtain 
; - V p= — Mol» (3.9) 
where we have used the condition that j is the gradient of a harmonic func- 
tion. Using eqs. (3.1), (3.2) and (3.9) it can be shown that the pressure, p, 


inside the fluid slab is given by 
p= Ey + Byz—} wipe + 1? cosh kz+ D, sinh kz 


_ Horio? (A, sinh kz+ A, cosh k=)| COS —, (3.10) 
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[S, sinh kd+ 8, cosh kd], 


Lae cosh 2kd 


where S,=a[(T+og+pdj2k?/K], | (3.12 a) 
; lt = 2 ‘2kk ky, By 4 
a Sa 5p [tei ea org 
where -T=surface tension of the fluid. 


(c) Dispersion relation _ 


Since we have assumed that the quantity a is periodic in time we can re- 
write the eq. (2.5) as 


—oaw* cos p=[—V p+ jxB—0g2Z]., (3.13) 


where the suffix zb on the right-hand side shows that the z-component of the 
quantities should be evaluated at the perturbed boundary. Using eqs. (3.2) to 
(3.8), (3.11) and (3.12) and including the effect of the homogeneous magnetic 
field produced by external sources, one can show, retaining terms of order a/d 
only, that eq. (3.13) reduces to ‘ 


' 


re Mojigkz(2kd—1) kk 4 
Ree Oo/0O™r Sa RG) / 
& Tah an ios tanh 2kd. Cea 


4. Discussion 


The dispersion relation given by eq. (3.14) is a generalisation of the one ob- 
tained by Rayleigh [6] to the case when the fluid is an electrical conductor 
carrying an uniform current in the direction of a homogeneous magnetic field 
produced by external sources. We see that if we put k,=0 in the eq. (3.14) 
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me a om Effi F 
re ve vector of the modes is < 
perpendicular to the direction of the current. 


- 


Sd = Oslom | Of tag oats ae 


ea 
the effect of both current and magnetic field is removed, but if we put k,=0 
n the eq. (3.14) the effect of the external magnetic field is removed and the 
effect of current is maximum because k,=k in this case. This is due to the fact 
that the fluid is assumed to have such low electrical conductivity that the 
magnetic field is not frozen inside. If B,=0 and we also neglect surface tension 
-and gravity eq. (3.14) shows that the fluid slab is stable for perturbations of 
wavelength 4<A, and unstable when A>A, where 


. 


Tie 


A,=4ad. (3.15) 


2 


a, wr 


Thus under these conditions the surface of the fluid is unstable when it carries 
-a volume current and deforms at a maximum rate into ridges and valleys 
directed along the current. When surface tension and gravity are also included, 
‘the critical wavelength, A,, at which the surface is neutral to perturbation, is 
larger than that given by eq. (3.15). 


f 


~ 


If we define a non-dimensional parameter 


a = 2Qkd, (3.16) 


‘the value of x=a, at which the surface is neutral to perturbation and the 
value x=2, with k,=0 at which it grows at a maximum rate when B,=0 are 
calculated graphically for mercury (7’=487 dyne em’; 9 =13.6 gm em~°*) for 
‘three values of current density and two values of depth and the results are 
shown in the Table 1. 

Corresponding to #=2,, the maximum growth rate of the unstable modes is 
calculated from eq. (3.14) when B,=0 and the results are shown in Fig. 2, 
‘where it will be seen that the characteristic time of growth of instabilities is 
of the order of a hundred ms and that this quantity is not very sensitive to 
the depth of the fluid tank. The values of current density used in Fig. 2 are 
of the order of a few hundred amp cm ™ which corresponds to a power of 
order 10* watts for a tank of mercury with depth 0.1 cm, length 100 cm and 


width 30 cm. 
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200 400 600 1000 
Acm 


Fig. 2. Variation of the growth rate of infinitesmal sinusoidal perturbations of the free sur - 

face of a mercury tank with electric current density. The wave vector of the perturbation 

is perpendicular to the direction of the current. Effects of surface tension and gravity are 

included. The value of the parameter 2d is equal to the depth of the tank, the bottom of 
which is assumed to be in contact with a rigid surface. 


If one puts 7, g and B, equal to zero in eq. (3.14) and considers the dis- 
turbance with k,=k, one can show that the e-folding time of the disturbance 
is approximately equal to 4d/V, where V, is the Alfvén wave velocity corre- 
sponding to the magnetic field at the boundary of the slab. This result is in 
conformity with the observation made by Murty (§ 4 (b) of [4]) that the 
time scales defined by Alfvén wave propagation over characteristic distances of 
the problem play a significant role even in finitely conducting media if the 
electric currents are maintained by external sources. 

If the value of B,+0, the behaviour of the fluid slab is more complicated. 
An examination of the dispersion relation shows that an external magnetic field 
directed along the current makes the fluid slab unstable when k, and k, have 
opposite signs. In such a case the wave vector of the perturbation lies in the 
second quadrant of the x-y-plane of the free surface and the resultant current 
density vector given by eq. (8.2) lies in the first quadrant of the x-y-plane 
at the points where the surface is perturbed upwards. This means that the wave 
vector of the unstable modes is so directed as to make the current density tend 
to be force free. This result is the same as the one obtained in the case of an 
infinite cylinder carrying axial current in the presence of a homogeneous magnetic 
field directed parallel to the current [4]. 

If the value of the external magnetic field, B,, is so large that the effect of 
all other terms in eq. (3.14) can be neglected except the one containing B,, 
we see that w° has numerically the largest value when |k,|=|k,|—k/ V2. There- 
fore, under these conditions the free surface is highly unstable for large wave- 
numbers of the perturbations with the wave vector lying in the second quadrant 
of the w-y-plane of the free surface inclined at an angle of 45° to the y-axis 
(direction of the current and the homogeneous magnetic field). 
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PART II 


= 


stability of a fluid slab with both the surfaces free, carrying a uni- 
3 form current in the presence of a homogeneous magnetic field 


G 


-_ 5. Formulation of the problem 

e 

_ In this section we shall consider the stability of an incompressible, inviscid 
and nonmagnetic fluid slab with both surfaces free and carrying a uniform cur- 
‘rent in the presence of an external homogeneous magnetic field. The geometry 
of this problem is the same as the one shown Fig. 1, except for the difference 
that we shall not include the influence of gravity and assume that the surface 
_2=—d is also a free surface. This configuration can approximately be realised 
by a jet of rectangular cross-section, the thickness of which is small compared 
with the dimensions of the jet in the other directions. As before we assume 
that there is a uniform current given by j=(0, 7), 0) flowing in the slab and 
that the current is closed symmetrically on either side of the slab at infinitely 
large distance. Also we introduce a uniform external magnetic field along the 
direction of the current. 

When the two surfaces of the fluid slab are free there are in general three 
possible ways in which a sinusoidal perturbation can deform these surfaces. The 
deformed surfaces can be (1) symmetric or (2) anti-symmetric or (3) be neither 
symmetric nor anti-symmetric with respect to the plane z=0. We shall show 
-that only the cases (1) and (2) arise when B,=0 and only the case (3) arises 
when B,+0. 

In order to do so we assume that the surface z= +d (hereafter called upper 

surface) and the surface z= —d (hereafter called lower surfase) are respectively 
deformed to shapes given by 


z=d-+a cos g, (5.1 a) 
and z= —d—aa cos 9, (5.1 b) 


where « is a nondimensional number. As before we assume that a@ is periodic 
in time (eq. (2.4)) and that a<d. We also demand that «a<d. The solution 
gives w and also the possible values of « which then determine the modes that 
-are most unstable. 

The basic equations that we use are again the set of eqs. (2.5) to (2.9) with 
eqs. (2.1) and (2.2) to describe the equilibrium conditions. 


> 


6. Solution for w 


The method of obtaining the value w in the present case is similar to the 
one adopted in §3 in the case of one free surface, except for differences arising 
out of the boundary conditions. The general form of the distribution of current 
density, magnetic field and pressure inside the fluid are given by eqs. (3.2), 
(3.7) and (3.10) respectively. The constants in eq. (3.2) have now to be deter- 
mined by demanding that there is no current density normal to the boundaries 
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Ta 
£=P+Q-M-N, 


H=O4N, (6.3) 
6-Q+M-P-N, (6.30) 


“where P= E T he? + 5 plielt ke (kd —cosh kd exp (— ka))| tanhkd, (64a 


= 57 i Bubs k, coth kd, (6.4b) 


M = Ja i+ spamoitlt (kd — sinh kdexp(—kd))| coth kd, (6.4¢ 


1 


and N= 57 ho Buksk, tanh kd. (6.4d) 


Eq. (6.2) shows that there are twe possible values of « and correspondingly 
there are two possible values of wm given by eq. (6.1), one for each value of a. 
The general form of the solution for @ is very involved as may be seen from 
eqs. (6.1) to (6.4). We therefore discuss the limiting cases when the external 
magnetic field is (1) zero, (2) very small and (3) very large, as these cases turn 
out to be particularly simple. 


7. Limiting cases 
(a) Bu=0 
If B,=0, we see from eqs. (6.4) that both Q and N are zero. Therefore n 
is also zero (eq. (6.3b)) and = —@ from eqs. (6.3). We then obtain from eq. 
(6.2) that under this limiting condition 
a=+1. (6.5) 
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In view of the general form of the perturbation we have started with (eq. (5.1)) 
we find that the fundamental modes are either symmetric («=1) or anti-sym- 


metric («= —1). The corresponding dispersion relations are, for symmetric per- 
turbation, : 
+ 3 = 9 
a GaP, (6.6) 
: Q 
and for anti-symmetric perturbation, 
2 
aw” =e M. (6.7) 


It can be shown that M20 for all values of kd while P can be negative for 
some range of values of kd. Therefore we conclude that the fluid slab carrying | 
a volume current in the absence of external magnetic field with both its sur- 
faces free can be unstable for symmetric perturbations only. Further we note 
that in the presence of electric current the dispersion relation for symmetric 
deformation (eqs. (6.6) and (6.4a)) is not identical to that for a fluid slab with 
half the thickness and with one surface rigid and the other remaining free (eq. 
(3.14) with g=0, B,=0 and replacing 2d by d). This is due to the fact that 
the two sides of the plane z=0 in the present case are not independent but are 
in interaction due to the magnetic field due to the current flowing in the fluid. 


(b) [Bu/MoJod] <1 


If the external homogeneous magnetic field is small compared with the field 
‘produced by the current at the boundary of the slab, it can be shown from 
eqs. (6.2) to (6.4) that the two possible values of «, «, and a, are given by 


er 


a,= +1+,+0("), (6.8 a) 

= —1+p,+ O(p'), (6.8 b) 

where B,= 6 coth kd, (6.9 a) 
and fB.=f tanh kd, (6.9 b) 
— 10 Bu ks fy 6.10 

where ares ar (6.10) 


Using eqs. ‘(6.b), (6.8) and (6.9) one can show that up to the terms of first 
order in B,/u)j,d the dispersion relations are the same as those given by eqs. 
(6.6) and (6.7). 


(c) [Bu/mojod]>1 
In the limit when the external homogeneous field is very large compared to 


the field at the boundary due to the current one can show from eqs. (6.2) to 


6.4) that 
. a= —exp (+2kd), (6.11) 
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a eq. (6.12) k, and ky may y bares cline le same 
Therefore gither of the values of « given by eq. (6.11) 
stable perturbations. The maximum instability occu 
tion when |k, |=|4|=4/V2 and the wave vector ‘gine 
of the x-y-plane, and the growth rate is proportional to (bE me 
- We see from eqs. (6.8) and (6.11) that «+1 only when B,+0, 
that the unstable modes of the fluid slab with both of its 
neither symmetric nor anti-symmetric with respect to she mid-y 
sence of an cea magnetic field. 
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